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ABSTRACT 

The measurement of the large-scale B-mode polarization in the cosmic microwave 
background (CMB) is a fundamental goal of future CMB experiments. However, be¬ 
cause of unprecedented sensitivity, future CMB experiments will be much more sen¬ 
sitive to any imperfect modelling of the Galactic foreground polarization in the re¬ 
construction of the primordial B-mode signal. We compare the sensitivity to B-modes 
of different concepts of CMB satellite missions (LiteBIRD, COrE, COrE-F, PRISM, 
EPIC, PIXIE) in the presence of Galactic foregrounds. In particular, we quantify the 
impact on the tensor-to-scalar parameter of incorrect foreground modelling in the com¬ 
ponent separation process. Using Bayesian htting and Gibbs sampling, we perform the 
separation of the CMB and Galactic foreground B-modes. The recovered CMB B-mode 
power spectrum is used to compute the likelihood distribution of the tensor-to-scalar 
ratio. We focus the analysis to the very large angular scales that can be probed only by 
CMB space missions, i.e. the Reionization bump, where primordial B-modes dominate 
over spurious B-modes induced by gravitational leasing. We find that fitting a single 
modihed blackbody component for thermal dust where the “real” sky consists of two 
dust components strongly bias the estimation of the tensor-to-scalar ratio by more 
than 5a for the most sensitive experiments. Neglecting in the parametric model the 
curvature of the synchrotron spectral index may bias the estimated tensor-to-scalar 
ratio by more than la. For sensitive CMB experiments, omitting in the foreground 
modelling a 1% polarized spinning dust component may induce a non-negligible bias 
in the estimated tensor-to-scalar ratio. 

Key words: cosmic microwave background - inflation - early universe - diffuse 
radiation - polarization - methods: analytical 


1 INTRODUCTION 

The detection of the CMB B-mode polarization, due to pri¬ 
mordial gravitational waves, would provide a definitive proof 
that inflation occurred in the Early Universe. This is chal¬ 
lenging because highly-polarized Galactic foregrounds dom¬ 
inate, by several orders of magnitude, the amplitude of the 
large-scale CMB B-modes. Therefore, the detection of the 
CMB B-modes is, by definition, a component separation 
problem (Betoule et al. 2009; Dunkley et al. 2009; Efstathiou 
et al. 2009; Bonaldi & Ricciardi 2011; Errard et al. 2011; Er- 
rard & Stompor 2012). 
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Ade et al. (2014) recently claimed the first detection 
of the primordial CMB B-modes at a tensor-to-scalar ra¬ 
tio r = 0.2 in the high-latitude region of the sky observed 
by the ground-based CMB experiment BICEP2. The fore¬ 
ground subtraction performed in the BICEP2 analysis was 
relying on poor assumptions on the polarized foreground 
model, which was constructed from the Planck intensity map 
of thermal dust (Planck Collaboration et al. 2014a). More 
recently, the Planck Collaboration published their first po¬ 
larized foreground data and showed that the BICEP2 de¬ 
tection of B-modes was likely to be due to the polarized 
thermal dust emission (Planck Collaboration et al. 2014b). 
In other words, the polarized dust signal in the BICEP2 field 
was in first instance significantly under-estimated. The joint 
BICEP2/Keck and Planck analysis has now shown this to 
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be the case with a combined upper limit of r < 0.11 for a 
pivot scale of ko = 0.002 Mpc“^ (BICEP2/Keck and Planck 
Collaborations et al. 2015). 

Due to the presence of the atmosphere, ground-based 
experiments, like BICEP2, are limited in frequency coverage. 
Furthermore, ground-based experiments are typically lim¬ 
ited in terms of sky area coverage. This may generate larger 
uncertainties to the B-mode signal detected by ground-based 
experiments, limiting their ability to achieve robust fore¬ 
ground subtraction on large angular scales (> 1° or multi¬ 
poles i < 200), where the cosmological signal from inflation 
is detectable. For example, the Reionization bump aut £ ~ 10, 
where the signal is due primarily to pure tensor perturba¬ 
tion modes. This is the reason why several fourth generation 
space-borne B-mode satellite concepts have emerged around 
the world (see Section 2). It is therefore highly likely that 
there will be at least one B-mode satellite in operation in the 
near future that will have the sensitivity to detect B-modes 
at the level of r ~ 10“® (e.g. Matsumura et al. 2013). 

Component separation is clearly one of the biggest chal¬ 
lenges for any B-mode satellite (Baccigalupi et al. 2004; 
Stivoli et al. 2006; Betoule et al. 2009; Dunkley et al. 2009; 
Bonaldi & Ricciardi 2011; Errard & Stompor 2012). Most 
of the component separation methods rely on prior physical 
assumptions on foregrounds to clean the CMB (see Planck 
Collaboration et al. (2014c) for a panel of methods used in 
the PZancfc temperature data analysis): Commander (Eriksen 
et al. 2008) performs a Bayesian fitting based on a paramet¬ 
ric model for foregrounds, SMICA (Delabrouille et al. 2003; 
Cardoso et al. 2008) does not rely on a parametric model 
for foregrounds but still assumes a given number of indepen¬ 
dent foreground components to perform the spectral match¬ 
ing between the data and the model. However, the Galac¬ 
tic foreground B-modes, unlike the CMB, have an unknown 
spectral signature and, even worse, the exact number of in¬ 
dependent polarized foregrounds is not known. More recent 
developments have been achieved in order to relax any prior 
assumption on the number of foregrounds, e.g. the GNILC 
method (Remazeilles et al. 2011). 

Given the unprecedented sensitivity of future B-mode 
satellites, any slightly incorrect foreground modelling in the 
component separation process is expected to strongly bias 
the reconstruction of the CMB B-mode signal (Armitage- 
Caplan et al. 2012). 

In this work we first compare the sensitivity of various 
CMB B-mode satellite concepts that have been put forward 
in recent years by evaluating the posterior probability dis¬ 
tribution of the key cosmological parameter, the tensor-to- 
scalar ratio r. We then go on to evaluate the impact of incor¬ 
rect foreground modelling on the estimation of r. Incorrect 
modelling of foregrounds investigated in this work includes: 
i) fitting for a single spectral component of thermal dust 
where the simulated sky data contain a second spectral com¬ 
ponent of thermal dust not in the data model, ii) neglecting 
the curvature of the synchrotron spectral index in the data 
model, and iii) omitting the spinning dust polarization in the 
data model. Our aim is to achieve a comprehensive study 
by implementing a Bayesian fitting component separation 
method. Commander, on a set of simulations with different 
foreground models and different instrument designs. Such 
a forecasting analysis is primordial for the optimisation of 
the concept of a future B-mode spatial mission. Note that 


in this comparative forecast analysis we do not consider the 
ability of each experiment to correct for lensing, limiting the 
analysis to low multipoles {£ < 12) where lensing has neg¬ 
ligible impact. In a realistic situation, space missions would 
also consider intermediate multipoles to guarantee a robust 
detection of the CMB B-modes, in which case lensing should 
be included for a fair comparison on larger range of multi¬ 
poles, and the ability of high-resolution experiments to cor¬ 
rect for lensing could make the difference. This should be 
investigated in a separate work. 

In Section 2 we summarise the main CMB B-mode satel¬ 
lites that have been proposed in recent times and tabu¬ 
late their instruments in terms of frequency channels, beam 
sizes and sensitivities in each channel. Section 3 discusses 
the sky models that are used in the simulations. Section 4 
summarises the method, based on the Commander Bayesian 
parametric fitting code. Results are given in Section 5 and 
conclusions summarised in Section 6. 


2 CMB B-MODE SATELLITE EXPERIMENT 
CONCEPTS 

A number of concepts for CMB B-mode satellite missions 
have emerged over recent years, particularly after the success 
of COBE and WMAP, and more recently, Planck (Planck 
Collaboration et al. 2011a). These proposed missions are fo¬ 
cussed on measuring the CMB polarization across a wide 
range of angular scales to an accuracy that is limited by 
astrophysical foregrounds. The sensitivities to the CMB B- 
modes are at least 2 orders of magnitude better than the 
current general of ground- and balloon-based experiments, 
which would allow a r ~ 0.001 to be detected. Each concept 
is briefly described below. Table 1 gives the basic experimen¬ 
tal set-up for each concept, based on the designs described 
in the literature. 

2.1 EPIC 

EPIC, Experimental Probe of Inflationary Cosmology (Bock 
et al. 2008, 2009), is a U.S. B-mode satellite mission concept 
proposed to the National Aeronautics and Space Admin¬ 
istration (NASA). The EPIC study team has investigated 
three possible designs with different instrumental character¬ 
istics. 

2.1.1 EPIC Low-Cost option with TES detectors 

The Low-Cost option (Bock et al. 2008), termed EPIC-LC, 
consists of six 30 cm refracting telescopes for an optimisa¬ 
tion of the cost of the mission, large arrays of Transition- 
Edge Superconducting (TES) bolometers, and a rotating 
half-wave plate in front of the optics to modulate the po¬ 
larization signal without modulating any polarization sys- 
tematics generated from the optics. It has a limited resolu¬ 
tion of between IGarcmin and 155arcmin and 7 frequency 
channels distributed between 30 GHz and 300 GHz. EPIC- 
LC is therefore focussed on the large-scale primordial CMB 
B-mode science. Two types of detectors were considered: i) 
Neutron-Transmutation-Doped (NTD) Ge Transition Edge 
Sensor (TES) bolometers with unmultiplexed JFET ampli¬ 
fiers (the same technology developed for the Planck HFI), 
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Concept name 

Leading country/ Frequencies 
institution [GHz] 

Beam size 

FWHM [arcmin] 

Sensitivities 
[fiK deg] 

Reference / notes 

EPIC-LC-TES 

U.S.A. 

(NASA) 

30,40,60,90 

135,200,300 

155,116,77,52, 

34,23,16 

0.460,0.156,0.085,0.037, 

0.035,0.037,0.062 

EPIC Low-Cost option with 

TES detectors (Bock et al. 2008) 

EPIC-CS 

U.S.A. 

(NASA) 

30,45,70,100, 

150,220,340,500 

15.5,10.3,6.6,4.6, 

3.1,2.1,1.4,0.9 

0.683,0.367,0.150,0.117 

0.117,0.183,0.883,7.50 

EPIC Comprehensive-Science option 
(Bock et al. 2008) 

EPIC-IM-4K 

U.S.A. 

(NASA) 

30,45,70,100,150 

220,340,500,850 

28,19,12,8.4,5.6 

3.8,2.5,1.7,1.0 

0.147,0.061,0.027,0.018,0.014, 

0.027,0.058,0.014,0.012 

EPIC Intermediate with 

4K mirror (Bock et al. 2009) 

LiteBIRD 

Japan 

(JAXA) 

60,78,100, 

140,195,280 

75,58,45, 

32,24,16 

0.172,0.108,0.078, 

0.062,0.0517,0.063 

(Matsumura et al. 2013) 

COrE 

Europe 

(ESA) 

45,75,105,135,165, 

195,225,255,285,315, 

375,435,555,675,795 

23.3,14,10,7.8,6.4, 

5.4,4.7,4.1,3.7,3.3, 

2.8,2.4,1.9,1.6,1.3 

0.150,0.078,0.077,0.075,0.077 

0.075,0.075,0.173,0.283,0.767, 

1.95,4.25,9.82,57.0,348.0 

ESA M mission concept 

(The COrE Collaboration et al. 2011) 

COrE + 

Light 

Europe 

(ESA) 

60,70,80,90,100, 

115,130,145,160,175, 

195,220,255,295,340, 

390,450,520,600 

21.0,18.0,15.8,14.0,12.6, 

11.0,9.7,8.7,7.9,7.2, 

6.5,5.7,5.0,4.3,3.7, 

3.2,2.8,2.4,2.1 

0.485,0.467,0.320,0.257,0.197, 

0.138,0.110,0.092,0.092,0.090, 

0.090,0.135,0.218,0.430,0.817, 

1.645,4.205,10.535,15.848 

ESA M mission concept 3^ 

COrE + 

Extended 

Europe 

(ESA) 

60,70,80,90,100, 

115,130,145,160,175, 

195,220,255,295,340, 

390,450,520,600,700,800 

14.0,12.0,10.5,9.3,8.4, 
7.3,6.5,5.8,5.3,4.8, 

4.3,3.8,3.3,2.9,2.5, 

2.2,1.9,1.6,1.4,1.2, 1.1 

0.342,0.233,0.160,0.123,0.098, 

0.073,0.057,0.057,0.057,0.058, 

0.063,0.090,0.152,0.220,0.422, 

0.790,1.982,5.632,20.05,93.5,203 

ESA M mission concept 4^ 

PRISM 

Europe 

(ESA) 

30,36,43,51,62, 

75,90,105,135,160, 

185,200,220,265,300, 

320,395,460,555,660 

17,14,12,10,8.2, 

6.8,5.7,4.8,3.8,3.2 

2.8,2.5,2.3,1.9,1.7, 
1.6,1.3,1.1,0.92,0.77 

0.211,0.141,0.133,0.103,0.098, 

0.093,0.078,0.068,0.061,0.0572 

0.059,0.061,0.064,0.073,0.085, 

0.092,0.135,0.197,0.404,0.953 

ESA L mission concept 
(Andre et al. 2014) 

PIXIE 

U.S.A. 

(NASA) 

30,60,90,120,150, 

180,210,240,270,300, 

330,360,390,420,450, 

480,510,540,570,600, 

630,660,690,720,750, 

780,810,840,870,900, 

930,960,990,1020,1050, 

1080,1110,1140,1170,1200 

96.0 (constant) 

5.180,1.390,0.691,0.454,0.352, 

0.307,0.292,0.297,0.319,0.358 

0.418,0.503,0.623,0.790,1.020, 

1.350,1.800,2.440,3.350,4.660, 

6.550,9.280,13.30,19.10,27.70, 

40.50,59.60,88.20,131.00,196.00, 

294.00,444.00,672.00,1020,1560, 

2390,3670,5660,8750,13600 

(Kogut et al. 2011) 


Table 1. Summary of proposed CMB B-mode satellite missions. The frequencies, beam sizes, and sensitivities are the nominal values 
for each of the missions (see text). The sensitivities are to Stokes Q (or U) and are given in CMB thermodynamic units. 


or ii) improved SQUID-multiplexed TES bolometers, which 
provide higher sensitivity. Here, we only provide results for 
the higher sensitivity TES option. 


2.1.2 EPIC Comprehensive-Science option 

Comprehensive-Science option (Bock et al. 2008), termed 
EPIC-CS, uses a 3m diameter dish, which guarantees ac¬ 
cess to a broader CMB science case compared to EPIC-LC. 
EPIC-CS has a higher resolution (« 5 arcmin in the CMB 
channel and up to 1 arcmin at high frequency) and a broader 
frequency range of observation (8 channels from 30 GHz to 
500 GHz) than the low-cost option. However, EPIC-GS has 
a lower sensitivity (0.12/rK deg at 150 GHz) than EPIC-LG 
(0.03/iK deg at 135 GHz). 


2.1.3 EPIC Intermediate option with 4K mirror 

The Intermediate option (Bock et al. 2009), termed EPIC- 
IM, consists of an intermediate size telescope (1.4 m diam¬ 
eter) cooled to 4K to reduce thermal noise. The angular 
resolution of EPIG-IM is similar to that of EPIG-CS, but 
the sensitivity of EPIG-IM is higher due to a mechanical 
4K cooling system and a focal plane of detectors designed 
in hexagonal cells. Moreover, EPIC-IM can observe the sky 
through a larger frequency range than the two other options 
with 9 channels distributed between 30 GHz and 850 GHz. 


There was also a less expensive option with a 30K mirror, 
which we do not consider here. 


2.2 LiteBIRD 

LiteBIRD, Lite (Light) satellite for the studies of B-mode 
polarization and Inflation from cosmic background Radia¬ 
tion Detection (Matsumura et al. 2013), is a Japanese satel¬ 
lite mission concept entirely dedicated to the detection of 
the large-scale primordial CMB B-mode polarization. This 
space mission is led by the Japan Aerospace Exploration 
Agency (JAXA). It was selected in March 2014 as one of 
27 highest-priority large projects by the Science Council of 
Japan. 

LiteBIRD is a light satellite with two small mirrors of 
« 60 cm diameter cooled to 4K, including 2000 TES (or 
MKIDs) bolometric detectors, and a half-wave plate in front 
of the instrument in order to modulate the polarized signal 
and to reduce the systematics. LiteBIRD should observe the 
sky within six frequency bands between 60 GHz and 280 GHz 
with a typical sensitivity of 0.06/rK deg at the CMB fre¬ 
quency channel (140GHz). 

Being focused on the detection of large-scale CMB B- 
modes, LiteBIRD does not rely on high-resolution chan¬ 
nel observations. The highest resolution is reached at the 
280 GHz channel with a 16 arcmin beam full width at half 
maximum (FWHM). 
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2.3 COrE 

COrE, Cosmic Origin Explorer (The COrE Collaboration 
et al. 2011) is a European B-mode satellite concept proposed 
to the European Space Agency (ESA) in 2011 as a medinm- 
class (M) mission. 

The design of COrE allows for observing the microwave 
sky in polarization over 15 frequency bands between 45 GHz 
and 795 GHz. While the COrE experiment provides an angu¬ 
lar resolution comparable to Planck, the sensitivity of COrE 
is 10 to 30 times better than Planck because of a hundred 
times more detectors. COrE includes a rotating half-wave 
plate in order to better control the systematics (in partic¬ 
ular temperature leakage) when measuring the polarization 
modes. Due to the presence of the half-wave plate, the size 
of the antenna can not be too large (« 1.5 m) in order to 
meet the size requirements of a M-class spatial mission. 

2.4 COrE+ 

COrE-b, Cosmic Origin Explorer Plus (see Table 1 for ref¬ 
erence) is a European B-mode satellite concept proposed to 
ESA in 2015 as an M-class mission. It was a follow-on from 
the original COrE concept. Two different designs have been 
investigated: a low-cost light satellite and an extended ver¬ 
sion with a larger telescope mirror. 

2.4-1 COrE+ Light concept 

In the COrE-b Light concept, the telescope has aim aper¬ 
ture allowing for Garcmin resolution at 200 GHz. COrE-b 
Light can observe the sky through 19 frequency channels 
between 30 GHz and 600 GHz with a claimed sensitivity of 
2.5/rKarcmin in CMB polarization after foreground clean¬ 
ing. COrE-b Light is the European low-cost (~ €550M) and 
low-risk option. 

2-4.2 COrE+ Extended concept 

The COrE-b Extended version has a 2 m aperture telescope, 
guaranteeing « 4arcmin resolution at 200 GHz, and an ex¬ 
tended frequency range of observations with 21 channels dis¬ 
tributed between 30 GHz and 800 GHz. The Extended ver¬ 
sion is expected to provide 1.5 pK arcmin sensitivity in GMB 
polarization after foreground cleaning. The Extended mis¬ 
sion was proposed to be funded with ESA and international 
partners, at a level of « €700M. 

2.5 PRISM 

PRISM, Polarized Radiation Imaging and Spectroscopy 
Mission (Andre et al. 2014), is a European space mission 
concept proposed to the European Space Agency (ESA) 
in 2013 as a large-class mission for characterizing the mi¬ 
crowave sky in intensity and polarization. 

With a relatively high resolution (few arcmin to 
1 arcmin), a high sensitivity (between 0.05/iK deg and 
0.9 pK deg), and a large frequency range through 20 observa¬ 
tion channels distributed from 30 GHz to 660 GHz, PRISM 
is designed to cover a broad range of GMB science. This 
includes GMB B-mode search, Sunyaev-Zeldovich cluster 
science, gravitational leasing. Cosmic Infrared Background, 


Galactic Science, and spectral distortions of the CMB abso¬ 
lute spectrum. PRISM is considered as the ultimate CMB 
mission in terms of sensitivity and resolution. 

PRISM consists in two instruments: a high angular res¬ 
olution polarimetric imager with a 3.5 m diameter telescope 
cooled to « 4 K, and a low angular resolution spectrometer 
to measure the spectral distortions to the absolute black- 
body spectrum of the CMB. For this study, we only simulate 
data from the imaging instrument. 

2.6 PIXIE 

The Primordial Inflation Explorer (PIXIE) is a U.S. space 
mission concept proposed to NASA, consisting of a polariz¬ 
ing Fourier transform spectrometer (FTS) with 400 spectral 
bands from 30 GHz to 6000 GHz (Kogut et al. 2011). Despite 
the fact that the instrument is designed to perform absolute 
spectroscopy, it can also be exploited to measure large-scale 
CMB B-mode polarization. The FTS operates as a nulling 
polarimeter by cancelling any unpolarized emission. PIXIE 
channels have approximately the same, relatively low, an¬ 
gular resolution of « 96 arcmin beam FWHM. The map 
sensitivity in each channel can be varied depending on the 
operation of the FTS, resulting in variable band widths. For 
this study, we use sensitivities appropriate for 40 channels 
between 30 and 1200GHz (A. Kogut, priv. comm.). 


3 SKY MODELS AND SIMULATIONS 

We consider a set of polarized sky simulations generated 
using version 1.7.8 of the Planck Sky Model (PSM) software 
package^ (Delabrouille et al. 2013). The software is also used 
to simulate Gaussian white noise at each frequency for each 
experiment, using the noise levels and angular resolutions 
given in Table 1. No bandpass effects are included, such 
that each frequency is assumed to be monochromatic. 

3.1 CMB 

The CMB is chosen to be a Gaussian field simulated from an 
input power spectrum defined by a set of standard cosmolog¬ 
ical parameters based {Qb = 0.0456, = 0.272, Qk = 0.0, 

Ho = 70.4kms“^Mpc"\ as = 0.809, A, = 2.441 x 10“®, 
Us = 0.963, T — 0.087) based on recent data®. We use three 
values of the tensor-to-scalar ratio, r, to test the effects of 
sample variance: 0.0 (no B-modes therefore no sample vari¬ 
ance), 0.001 (minimal sample variance, close to the sensitiv¬ 
ity level of the experiments we are testing) and 0.05 (just 
below the current upper limits of r < 0.1; BICEP2/Keck 
and Planck Gollaborations et al. 2015). Note that we use 
a pivot scale of ko = 0.05 Mpc“^ but the results are very 
weakly dependent on this because we have no tilt in the 
tensor spectrum (nr = 0.0). The input theoretical power 

® Note that during the last stages of the completion of this work, 
new results from the observations of the Planck satellite have 
been posted, pointing to a lower value of the optical depth to 
Reionization with r = 0.066 (Planck Collaboration et al. 2015e). 
This means that for a given tensor-to-scalar ratio, the amplitude 
of CMB B-modes with respect to foreground B-modes should be 
even lower. 
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spectrum is then calculated using the CAMB software (Lewis 
et al. 2000). We do not include the effects of gravitational 
lensing, which would contribute significant power at scales 
I > 100. For a fair comparison, we use the same realiza¬ 
tion of the CMB sky for each experiment and model. The 
main parameter that we vary is the tensor-to-scalar ratio, 
r. The CMB fluctnations in both intensity and polarization 
are scaled with freqnency according to the same spectral 
emission law, a{v), given by the derivative with respect to 
temperature of a blackbody spectrum at Tcmb = 2.725 K: 

( 1 ) 

'r=TcMB 


l{v) = 


dB,{T) 

dT 


3.2 Foregrounds 

The simulation of the polarized sky includes the emission 
from Galactic foregrounds: synchrotron radiation and ther¬ 
mal dust radiation. For specific purposes in the present work, 
we also include in some simulations of the set the polariza¬ 
tion radiation from a spinning dust component. We do not 
include the effects of polarized sources, since they are only 
dominant on smaller angular scales (£ > 100). 


3.2.1 Synchrotron 

For the synchrotron component, the and po¬ 

larization maps are simulated at different frequencies by ex¬ 
trapolating the WMAP polarization observations at 23 GHz 
(Kogut et al. 2007) through a simple power-law 

QT^ = Q2^GH.[^y/ 

= U23CH.{Q^/ ( 2 ) 

considering either a uniform spectral index, (3s = —3, or a 
non-uniform synchrotron spectral index (Miville-Deschenes 
et al. 2008) over the sky. We also consider for some of the 
simulations the possibility that the synchrotron spectral in¬ 
dex includes a frequency dependence, i.e. 

/3.(^) = -3 + C71og(^), (3) 

through the addition of a curvature term, C. We adopt val¬ 
ues of C = 0.3 and C = —0.3, resulting in a flattening 
or steepening synchrotron spectrum at frequencies above 
23 GHz (Kogut et al. 2007). 


3.2.2 Thermal dust 

For the thermal dust, the and polarization maps 

are derived from the intensity map as 

Qdust ^ cos {2'yd), 

jjdust ^ ^ 2 ^^) . ( 4 ) 

The dust polarization angle map, jd, and the geometric de¬ 
polarization map, Qd, due to magnetic field configuration 
(Miville-Deschenes et al. 2008), are coherent with those of 
the polarized synchrotron model. The dust polarization fac¬ 
tor is set to fd = 0.15, which, after correction by the geo¬ 
metric depolarization factor gd, results in an observed po¬ 
larization fraction, fdQd, of about 5% on average over the 


sky. This is slightly lower than the typical value at high lati¬ 
tudes (« 8 %) that has recently been observed by the Planck 
satellite at 353GHz (Planck Collaboration et al. 2015c, b); 
the precise values depend on the zero levels in the Planck 
data and the polarization fraction varies significantly over 
the sky up to a maximum of ~ 20 % along some sight-lines. 

For the dust intensity map, 1^, we choose to extrapolate 
at different frequencies either the Planck thermal dust opac¬ 
ity map at 353 GHz (Planck Collaboration et al. 2014a) or 
the Schlegel-Finkbeiner-Davis (SFD98) dust intensity map 
at 100 micron (Schlegel et al. 1998). We consider two al¬ 
ternative models of emission law for extrapolating the dust 
template over frequencies: either a single modified blackbody 
(MBB) spectrum (Model la and Model 2a in Table 2) 

Iiy = I 353 Bv{Td), (5) 

with uniform dust spectral index = 1.6 and uniform dust 
temperature Td = 18 K, in which case we extrapolate the 
Planck dust model at 353 GHz, I 353 , or a more complex 
model of emission (Model lb and Model 2b in Table 2) 
parametrized by two modified blackbodies based on model 
7 of Finkbeiner et al. (1999): 


P=I 


100/J.m 




accounting for different populations of dust grains: a fraction 
fi of cold dust component with temperature Ti = 9.6 K and 
spectral index (3i = 1.5, and a dominant fraction /2 of hot 
dust component with temperature T 2 = 16.4 K and spec¬ 
tral index (32 = 2.6. In this case, the SFD98 dust intensity 
map at 100 micron, Iiooum, is extrapolated over frequen¬ 
cies. We note that recent results from the Planck satellite 
suggest that the Rayleigh-Jeans thermal dust spectral index 
varies relatively little over the sky, with an average value of 
(3d ~ 1.59 ±0.02 and a la dispersion of 0.16 (Planck Collab¬ 
oration et al. 2015f). Therefore our model is arguably more 
complex than the real sky. However, this is partly compen¬ 
sated by the fact that we are using a smaller average polar¬ 
ization fraction than observed. These details do not change 
the major conclusions of this paper. 


3.2.3 Spinning dust 

For simulations Model If and Model 2f (Table 2), we also in¬ 
clude a spinning dust component with 1% fraction of polar¬ 
ization (Dickinson et al. 2011; Rubifio-Martfn et al. 2012a) 
and same polarization angles than thermal dust. The in¬ 
tensity map for the spinning dust component is taken from 
the Planck thermal dust map at 353 GHz (Planck Collabo¬ 
ration et al. 2014a), rescaled by 0.91 K/K using the corre¬ 
lation coefficient from Planck Collaboration et al. (2015d), 
and extrapolated at each frequency from the 23 GHz value 
assuming a cold neutral medium (CNM) model of emission 
law (Ali-Haimoud et al. 2009). 


^ http://conservancy.umn.edu/handle/11299/169642 
^ http: //www. ape .univ-paris7 . f r/-delabrou/PSM/psni.html 
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Model ID 

CMB 

Synchrotron 

Thermal dust 

Spinning dust 

Model Ob 

r ^ 0.05 

N/A 

N/A 

N/A 

Model la 

r ^ 0.05 

Uniform 

One MBB 

N/A 

Model lb 

r ^ 0.05 

Uniform 

Two MBB 

N/A 

Model Ic 

r ^ 0.05 

Non-uniform 

One MBB 

N/A 

Model Id 

r ^ 0.05 

-|-0.3 curvature 

One MBB 

N/A 

Model le 

r - 0.05 

— 0.3 curvature 

One MBB 

N/A 

Model If 

r - 0.05 

Uniform 

One MBB 

1% 

Model Oc 

r = 0.001 

N/A 

N/A 

N/A 

Model Oa 

r = 0 

N/A 

N/A 

N/A 

Model 2a 

r = 0 

Uniform 

One MBB 

N/A 

Model 2b 

r = 0 

Uniform 

Two MBB 

N/A 

Model 2c 

r = 0 

Non-uniform 

One MBB 

N/A 

Model 2d 

r = 0 

-|-0.3 curvature 

One MBB 

N/A 

Model 2e 

r = 0 

— 0.3 curvature 

One MBB 

N/A 

Model 2f 

r = 0 

Uniform 

One MBB 

1% 


Table 2. Summary of sky model simulations used in this work. 


3.3 Set of simulations 

The set of sky models is summarized in Table 2. For example, 
Model lb includes CMB B-modes with tensor-to-scalar ratio 
r = 0.05, and is a nominal foreground polarization model: 
basically the thermal dust is parametrized by two modified 
blackbodies (Finkbeiner et al. 1999), and the synchrotron 
by a uniform power law. 

The nominal noise levels in each frequency channel for 
the variety of B-mode satellite concepts are listed in Table 1 
and co-added to the sky simulations. We assume that the 
noise is thermal (no 1/f noise) and uniformly distributed 
on the sky. 

The alternative descriptions of the Galactic foreground 
polarization proposed in this set of simulations allows us to 
study the impact on the CMB B-modes of fitting an incor¬ 
rect model for thermal dust, synchrotron, or spinning dust. 
Specifically, we test the impact of under-estimating the com¬ 
plexity of the thermal dust frequency dependence by fitting 
to the data a single MBB dust component whereas the sky 
consists of two MBB dust components, of neglecting the 
synchrotron curvature or the variability of the synchrotron 
spectral index, and of omitting low-polarized Galactic fore¬ 
grounds such as the spinning dust radiation. 


4 METHOD 

4.1 Separation of the component maps 


zz, is a linear combination of the CMB component, the ther¬ 
mal dust emission and the synchrotron emission: 


d(j>, v) = 


a{u)s^^\p) 


+ —T 


\P) 


+ 

+ 


1 / bAtRp)) 


(^) 




(Td(p)) 


dust / \ 

s (p) 


( 7 ) 


We assume a uniform Gaussian instrumental noise n(p, u) 
for each experiment. Therefore, the fitting model is 
parametrized by 

s = = (/5,,Td,/3.) ,C7, = (8) 

where s collects the Q, and U maps of the CMB the 

thermal dust and the synchrotron The Galactic 

foregrounds are also parametrized by the spectral indices, 
Pd for the thermal dust and Ps for the synchrotron, and 
the temperature Tz for the thermal dust. Ce is the angular 
power spectrum of the CMB maps and can be written as 



/ cr 

r^TE 

cr \ 


Cl = 1 

1 

1 

/^EE 

cr ■ 

( 9 ) 

1 

\ cr 

cf^ 

cr ) 



While the CMB and the noise can be accurately 
parametrised from known theoretical models. Galactic fore¬ 
ground parametrisation can be subject to incorrect model¬ 
ling due to imperfect knowledge of the foreground physics. 
In the present work, we play with the second and third lines 
of Eq. 7 allowing for some mismatch between the fitted fore¬ 
ground model and the foreground data from the sky simu¬ 
lations for different B-mode experiments. 

Using Bayes’ theorem, the joint CMB-foreground pos¬ 
terior distribution is given by: 

P{s,f3,Ce\d) cxP{d\s,l3,Ce)Pis,l3,Ct) (10) 

where P (s, f3, Ce) is the prior distribution. The Commander 
algorithm performs a Gibbs sampling on the Monte-Carlo 


Markov chains (MCMC) as 




p[s\cf,r\d) , 


^(i + 1) 







(11) 


which converges to the sampling from the joint posterior 
P{s,/3,Ce\d). 

From the marginalized distribution of the CMB power 
spectrum 

P{Ci\d) = J P {s,/3,Ce\d) ds^'^^ds‘^^‘‘^ds‘«"^dPddPs, (12) 

the mean CMB power spectrum 


We perform the component separation and the reconstruc¬ 
tion of the CMB map and power spectrum by implement¬ 
ing a Bayesian parametric fitting method using Gibbs sam¬ 
pling through the implementation of the Commander algo¬ 
rithm (Eriksen et al. 2008). 

The model of the sky d(p,u), that we fit to the simu¬ 
lated data in direction p (pixel) in the sky and at frequency 


{Ce) = J P{Ci\d)CedCi 

therefore reduces to a simple ensemble average over the 
Gibbs samples: 

{Ci) (13) 

where N is the number of samples in the Markov chain. 
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Similarly, the covariance matrix of the Ce can be computed 
as the ensemble average over the Gibbs samples as 


1 ^ 2 


(14) 


Since the sky is assumed Gaussian and isotropic, the 
conditional distribution used in Eq. 11 to sample Ct is an 
Inverse-Gamma distribution when interpreted as a function 
of G,: 

(2^ + 1) / 1 I 

p 2 Ci (21+1 

P [Ct I s) oc - 




C\ 


{ 2 l +•^)/2 


(15) 


for which a simple textbook algorithm exists 
(Eriksen et al. 2004). The conditional distribution used to 
sample the amplitude of the sky components is a Gaussian 
distribution 


P{s\Ct,d) oc P{d\s,Ci)P{s\Ci) 

d-l/2)(d-s)^JV 


OC 

oc 




(-l/2)(s-s)^(S '^+N ^)(s-s) 


(16) 


where S and N are the covariance matrices in pixel space 
of the GMB and the noise respectively, and s is given by the 
Wiener filtered data 


= {S~^ +N~/ N~^d. 


(17) 


The sample can therefore be derived from the sample 

as the conjugate gradient solution of 


{S~'^ + N-/ s = N- 


(18) 


where uiq, Wi ~ A/’(0,1). 


4.2 Likelihood estimation of the cosmological 
parameters 

Once CMB power spectra and have been esti¬ 

mated from Eq. 13, as well as the covariance matrix from 
Eq. 14, we can perform a likelihood analysis of the cosmo¬ 
logical parameters r (optical depth to Reionization) and r 
(tensor-to-scalar ratio). 

At low multipoles, e.g. £ < 12, we make the following 
assumption on the CMB power spectrum: 


< 12 oc r , 


ryBB 

Cl < 12 OC r. 


(19) 


We then compute the likelihood distribution 

- 21n£ (Gr|r,r) oc [Oi - C/ (T,r)] E”' [g^ - C/ (r,r)] (20) 

for different values of r and r, i.e. for different ampli¬ 
tudes of a theoretical template of the CMB power spec¬ 
trum, Cf^ {T,r), (Eq 19). This allows us for drawing the 
likelihood distributions P{t) and P{r) respectively for the 
optical depth, r, and the tensor-to-scalar parameter, r. 

While r can be constrained by the amplitude of E- 
modes, the tensor-to-scalar ratio r relies on the sole ampli¬ 
tude of the CMB B-mode power spectrum, which unlike the 
amplitude of E-modes is many orders of magnitude smaller 
than the amplitude of the Galactic foreground polarization 
signal. For this reason r is much more sensitive than r to any 
imperfect foreground modelling. In the following, we show 
the impact of foregrounds on the tensor-to-scalar parameter 
only. 


5 RESULTS 

In this section we compare the sensitivity to large-scale B- 
modes (tensor-to-scalar ratio) of different satellite concepts, 
GOrE, COrE-b, LiteBIRD, PIXIE, EPIC, and PRISM, first 
assuming no Galactic foreground, then considering Galac¬ 
tic foreground cleaning with the parametric fitting method 
described in Section 4.1. In particular, we show the impact 
on the estimated tensor-to-scalar ratio of incorrect assump¬ 
tions on the Galactic foreground model. For all GMB ex¬ 
periments, we compute both the value from parametric 
fitting and the bias on the recovered tensor-to-scalar ratio 
resulting from incorrect foreground modelling. 

We consider different mismatches between the fore¬ 
ground model and the simulated data: i) incorrect spectral 
modelling of the thermal dust through the assumption of a 
single MBB spectrum instead of two MBB spectra, ii) im¬ 
pact of high-frequency channels and Galactic masking, iii) 
incorrect spectral modelling of the synchrotron through the 
omission of a curvature term in the spectral distribution or 
through the incorrect assumption of a uniform spectral in¬ 
dex, and iv) omission of an extra polarized component like 
spinning dust. 

The different sky simulations for polarized foregrounds, 
investigated in this work, are listed in Table 2. Regardless 
of the different sky simulations, we intentionally fit with 
Commander a unique parametric foreground model to the 
data given by Eq 7: a single MBB component for the thermal 
dust, and a uniform power-law with zero curvature for the 
synchrotron. We adopt Gaussian priors for the dust spectral 
index Pd = 1.6±0.3, the dust temperature Td = 18±0.05K, 
and the synchrotron spectral index Pa = —3±0.1. The Gaus¬ 
sian priors are multiplied by the Jeffreys prior to take into 
account the volume of likelihood space for non-linear param¬ 
eters (Eriksen et al. 2008). Therefore, we are considering a 
possible mismatch between the model and the data regard¬ 
ing the spectral behaviour of the thermal dust and the syn¬ 
chrotron, and also regarding the omission of the spinning 
dust polarization. 

The likelihood distributions of the tensor-to-scalar ratio 
for the different CMB satellite concepts and the different sky 
configurations considered in this work are plotted in Fig. 1 
(dust modelling), in Fig. 7 (synchrotron modelling), and in 
Fig. 8 (spinning dust modelling). The maximum likelihood 
value of r obtained for the different CMB satellite concepts 
and the different sky configurations, taken together with the 
X^-values testing the quality of the Bayesian parametric fit¬ 
ting at the component separation stage, are listed in Ta¬ 
bles 3, 5, and 6, along with their Icr errors. 


5.1 Overall sensitivity on large scales 

5.1.1 No foregrounds, r = 0 

In the absence of any CMB B-modes (r = 0) and any Galac¬ 
tic foreground (model 14 in Table 2) the simulated observa¬ 
tions of the sky at any frequency for any CMB experiment 
would be characterized by the combination of pure CMB 
E-mode polarization and instrumental noise. 

Such a simple configuration can provide forecasts on the 
instrumental sensitivity to B-modes of the different CMB 
satellite concepts considered in this work. In this case, the 
implementation of the Commander algorithm reduces to a 
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simple denoising of the CMB polarization in each frequency 
channel. The complete set of frequency channels is used in 
the Bayesian fitting: 15 channels for COrE from 45 GHz 
to 795 GHz, 19 channels for COrE+ Light from 60 GHz to 
600 GHz, 21 channels for GOrE+ Extended from 60 GHz to 
800 GHz, 6 channels for LiteBIRD from 30 GHz to 280 GHz, 
33 channels for PIXIE from 30 GHz to 990 GHz (we do 
not consider channels beyond 1000 GHz that are available 
to PIXIE since this would involve more complicated fore¬ 
ground modelling), 7 channels for EPIG-CS from 30 GHz to 
300 GHz, 8 channels for EPIC-GS from 30 GHz to 500 GHz, 9 
channels for EPIC-IM-4K from 30 GHz to 850 GHz, 20 chan¬ 
nels for PRISM from 30 GHz to 660 GHz. Depending on the 
different frequency coverages and the different instrumental 
sensitivities in each individual channel, the resulting over¬ 
all sensitivity to B-modes may vary from one CMB satellite 
concept to another. 

The r.m.s. of the B-mode polarization map resulting 
from pure instrumental noise residuals after component sep¬ 
aration provides the overall sensitivity to B-modes for each 
CMB experiment. This exercise is interesting as a pure 
test of sensitivity because in the absence of CMB B-modes 
the error is no longer dominated by cosmic variance. The 
overall sensitivity to B-modes for all CMB satellite con¬ 
cepts is shown in the top left panel of Fig. 1. The exper¬ 
iments PRISM, EPIC-IM, and COrE-|- Extended show the 
best instrumental sensitivity to B-modes while PIXIE and 
EPIC-CS show the largest r.m.s. noise in terms of B-mode 
polarization. It is interesting to note that, in the absence of 
any foreground and any lensing contamination, COrE and 
LiteBIRD have similar sensitivity to B-modes. This of course 
must be reconsidered in presence of foreground contamina¬ 
tion (see hereafter) and/or lensing contamination; having 
more frequency channels than LiteBIRD, COrE is better 
designed than LiteBIRD to reduce uncertainties on r af¬ 
ter foreground cleaning. The reduced number of frequency 
channels from LiteBIRD may reduce its ability to clean con¬ 
tamination by non-trivial foregrounds. In addition, high- 
resolution experiments, like COrE, should be in a better 
position than low-resolution experiments, like LiteBIRD, to 
correct for lensing B-mode contamination and reduce the 
uncertainty on r due to lensing. 


5.1.2 No foregrounds, r = 0.05 

When the sky includes CMB B-mode polarization with a 
tensor-to-scalar ratio r = 0.05 but still no Galactic fore¬ 
ground - Model Ob in Table 2 - we do not need to fit 
for Galactic foregrounds at the component separation stage 
with Commander but only for CMB and noise. The complete 
set of frequency channels of the CMB satellite concepts is 
used in the Bayesian fitting. The low-multipole CMB B- 
mode power spectrum resulting from component separation 
can be used afterwards to compute the likelihood distribu¬ 
tion of the tensor-to-scalar ratio (Eqs 19 and 20) by using 
multipole modes up to the Reionization bump {i < 12). 

As expected, in the absence of any Galactic foreground 
and a fortiori the absence of any assumptions about fore¬ 
grounds, we get an unbiased estimate of the recovered 


tensor-to-scalar ratio (top right panel of Fig. 1)."* In the 
absence of foregrounds, a r = 0.05 signal is detected at 8cr 
by PIXIE and at 9cr by all the other experiments. When 
the sky includes CMB B-mode polarization with a tensor- 
to-scalar ratio r = 0.05 but still no Galactic foregrounds, the 
error on the recovered tensor-to-scalar ratio is dominated by 
cosmic variance instead of instrumental noise for all CMB 
satellite concepts considered here. This is illustrated in the 
top right panel of Fig. 1 where the likelihood distribution 
of the tensor-to-scalar ratio P{r) show equal width for all 
CMB satellite concepts, despite the different instrumental 
characteristics of each satellite (also see Table 3 where the 
la errors on r have similar values for all CMB experiments). 
The sample variance can be reduced by analysing to higher 
Gvalues, however, this would then require de-lensing and can 
be achieved from ground- or balloon-based experiments. 


5.1.3 No foregrounds, r = 0.001 

Because of the high level of sensitivity of the CMB experi¬ 
ments considered here, even if we include in the likelihood 
analysis multipoles up to .^ < 50, the error on the estimated 
tensor-to-scalar ratio continues to be dominated by cosmic 
variance in the absence of any foreground, at least when 
r = 0.05. 

In the case where the amplitude of the CMB B-modes 
would be smaller on large scales, say r = 0.001, and there is 
still no foreground - Model Oc - the proportion in the error 
budget due to cosmic variance becomes mitigated with re¬ 
spect to the intrinsic noise of the experiments. In this case, 
each CMB experiment show distinct Icr errors on the recov¬ 
ered tensor-to-scalar ratio (Fig. 2 and Table 4). In particular, 
the experiments PRISM, COrE-l- Extended, and EPIC-IM, 
yield the smallest error on the tensor-to-scalar ratio estimate 
with a detection of a r = 0.001 signal at 6cr in the absence 
of foreground/lensing, therefore confirming the sensitivity 
results presented in the top left panel of Fig. 1. 

These numbers should be considered with caution be¬ 
cause in this forecast analysis we neglect any additional error 
caused by lensing B-modes. Especially when r is low, the 
impact of lensing might dominate the error on r over the 
contribution from both instrumental noise and cosmic vari¬ 
ance. Depending on the resolution of each experiment, the 
sensitivity to r after delensing may vary from an experiment 
to another (Smith et al. 2012). 


5 . 1.4 Galactic foreground B-modes 

The unprecedented sensitivity of next-generation CMB 
satellite missions is required to detect the primordial CMB 
B-mode signal on large scales, in particular if the energy 


^ Note that the template power spectrum used in the likelihood 
Eq. 20 is not the CAMS power spectrum with the true value r = 0.05 
but the power spectrum of the input CMB map realization of 
the simulation, for which the actual value of the realization is 
r = 0.042. In all the plots of the paper, we deliberately rescale 
the realization value r = 0.042 to the true value r = 0.05 in order 
to highlight biases resulting from incorrect foregound modelling 
not from a particular CMB realization. 
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Figure 1. Recovered posterior distribution P{r) of the tensor-to-scalar ratio and impact of incorrect dust modelling. The theoretical 
input tensor-to-scalar value (vertical solid black line) is r = 0 in the left-hand panels and r = 0.05 in the right-hand panels. Top panels: 
no foregrounds {left: Model Oa, right: Model Ob). Middle panels: correct foreground modelling {left: Model 2a, right: Model la). Bottom 
panels: incorrect spectral modelling of thermal dust {left: Model 2b, right: Model lb). Recovered tensor-to-scalar distributions: COrE 
(solid yellow), COrE-|- Light (solid light-blue), COrE-|- Extended (solid blue), LiteBIRD (dotted red), PIXIE (dashed green), EPIC-LC- 
TES (long-dashed yellow), EPIC-CS (long-dashed purple), EPIC-IM-4K (long-dashed orange), PRISM (dash three-dot black). The top 
left panel compares the overall sensitivity of the different satellites in the absence of foregrounds by showing for Model Oa (r = 0, no 
foregrounds) the r.m.s of the residual noise B-mode map after component separation. 


scale of Inflation is small, resulting in r « 0.1. As a con¬ 
sequence, the uncertainty in the estimation of the tensor- 
to-scalar ratio will no longer be dominated by instrumental 
noise on large scales. Rather, the main uncertainty on r will 
come from the contamination by highly-polarized Galactic 
foregrounds. 

Considering the presence of Galactic foreground polar¬ 
ization in the observations, we can successfully perform fore¬ 


ground cleaning / component separation by implementing 
the Bayesian parametric fitting method described in Sec¬ 
tion 4.1, as long as our fitting model correctly matches the 
real (simulated) data in terms of parametrization of Galactic 
foregrounds. 

Model la includes Galactic foreground polarization 
through synchrotron radiation with a uniform power-law 
spectral distribution {fis = —3) and thermal dust radia- 
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tion characterized by a single MBB emissivity {fid = 1-6, 
Td = 18K), in addition to CMB polarization with r = 0.05 
and instrumental noise. In this case, we perform compo¬ 
nent separation by fitting a parametric model that perfectly 
matches the simulated sky, i.e. single MBB component for 
thermal dust and synchrotron power-law with uniform spec¬ 
tral index over the sky (Gaussian priors Ps ~ —3 ± 0.1, 
Pd = 1.6 ± 0.3, Td = 18 ± 0.05K). Moreover, the complete 
set of frequency channels of each experiment is used for the 
fit. Model 2a is identical except that there are no CMB B- 
modes (r = 0). 

The panels on the second row of Fig. 1 show that, with 
no modelling error on the Galactic foreground parametriza- 
tion (Models la and 2a), the tensor-to-scalar ratio is suc¬ 
cessfully recovered without any bias for all the CMB satel¬ 
lite experiments. By comparing the right panels of Fig. 1 
from the first row (Model Ob, no foreground) and the sec¬ 
ond row (Model la, including foregrounds), we can see that 
the presence of Galactic foregrounds inflates the error of the 
recovered tensor-to-scalar ratio. From the first two blocks of 
Table 3 (Models Ob and la) we can measure that the error 
inflation due to foreground contamination goes from « 28 % 
(i.e. the error is increased by a factor 1.28) for PRISM to 
« 40 % (i.e. the error is increased by a factor 1.4) for PIXIE. 
In presence of foregrounds, we find after component sepa¬ 
ration that r = 0.05 is now detected at 6 a for PIXIE and 
EPIC-CS and detected at 6.5a by LiteBIRD, COrE-l- Light 
and EPIC-LC, while it is detected at 7a by COrE, COrE-f 
Extended, EPIC-IM, and PRISM. 

The left panel on the second row of Fig. 1 (Model 2a, 
r = 0) shows different widths for the likelihood distribu¬ 
tion of the tensor-to-scalar ratio from each CMB experiment. 
This illustrates the ability of each experiment in removing 
the polarized foregrounds from the observations. In particu¬ 
lar, PIXIE shows less sensitivity than the other experiments 
because it is intrinsically less sensitive than the other ex¬ 
periments as shown in Sect. 5.1.1, with a la error on r of 
(Jr = 4.2 X 10“"^. Despite the large number of high-frequency 
channels of PIXIE, the uncertainty on r is still surprisingly 
larger than for the other experiments. Even if we include 
in the fit the whole set of frequency channels of PIXIE (> 
1 THz), yet the uncertainty on r is not significantly reduced. 
Actually, we find that even if the thermal dust is best fitted 
with PIXIE channels, nevertheless in this particular config¬ 
uration and simulation of the sky it is not the dust but the 
synchrotron contamination that dominates the increase in 
the uncertainty on r. Therefore, extra low-frequency chan¬ 
nels would better help than high-frequency channels in re¬ 
ducing the uncertainty on r. In particular, by combining 
the sensitive low-frequency C-BASS data at 5 GHz (Irfan 
et al. 2015) with PIXIE data we estimate that the error on 
r = 0.05 would be reduced by a factor 1.4, therefore com¬ 
pensating the increase of error due to synchrotron contami¬ 
nation. GOrE-h Extended shows the best overall sensitivity 
and ability in removing the foreground contamination with 
(Jr = 1.2 X 10“"*. The recovered values for r and their uncer¬ 
tainties are given in the fifth block of Table 3. 

In principle, we can measure the quality of the Bayesian 
fit / component separation with the distribution: 


X^{P) = 

U 


/ d{p, v) - m{p, u) 

V 


2 


( 21 ) 


where the sum runs over frequency channels, (Ji/(p) rep¬ 
resents noise in the data, and {d(p,v) — m(p,v)) is the 
mismatch between the observation data, d(p,v), and the 
parametric model, m(p,v), based on the spectral parame¬ 
ters. For a perfect fit, this quantity should be proportional 
to the number of degrees of freedom used in the fit, i.e. 
d.o.f = {Nch — Npar) where Nch is the number of observa¬ 
tion channels (twice the number of frequencies because of Q 
and U polarization observations) and Npar the number of 
parameters of the model. However, the effective number of 
parameters is not clearly defined in this Bayesian framework 
because of the presence of priors constraining the parameter 
space (Planck Gollaboration et al. 2014c). For this reason, 
it is difficult to compute the exact reduced chi-square dis¬ 
tribution, y^/d.o./, for the goodness-of-fit. However, we can 
use in this work the y^-values of Model la and Model 2a 
(Table 3) as the baseline for Hi® perfect fit of the fore¬ 
grounds. Similarly, the tensor-to-scalar ratio values of Model 
la and Model 2a (Table 3) can be regarded as the baseline 
values indicating correct foreground modelling. 

We can see that the mean x^ of the fit in the presence 
of Galactic foregrounds (Model la in Table 3) have simi¬ 
lar values that in the case with no foregrounds (Model Ob 
in Table 3). Any departure from the baseline y^-values of 
Model la and Model 2a indicates incorrect foreground mo¬ 
delling, which may result in a bias in the recovered tensor- 
to-scalar ratio. Gonversely, if the recovered tensor-to-scalar 
ratio shows a bias but the x^-Mslne remains close to the 
baseline value, then the fitting is still good over the range of 
frequencies considered and the bias on the tensor-to-scalar 
ratio is due to a lack of some frequency channels for the 
experiment. In this case, more frequency channels / data 
points would be needed for fitting a more complex fore¬ 
ground model, taking into account more degrees of free¬ 
dom (e.g. non-uniform spectral index, synchrotron curva¬ 
ture, etc.). 


5.2 Impact of incorrect foreground modelling 

As a counterpart of high sensitivity, next-generation GMB 
satellite experiments will be more sensitive to any incor¬ 
rect assumptions on the Galactic foreground parameters in 
any component separation algorithm. In particular, incor¬ 
rect foreground modelling may dramatically bias the estima¬ 
tion of tensor-to-scalar ratio. Thermal dust modelling might 
be critical in the high-frequency regime (> 100 GHz) of the 
observations, while synchrotron modelling and spinning dust 
polarization may play a non-negligible role at low frequency 
(< 100 GHz) for the accurate estimation of the tensor-to- 
scalar ratio. 


5.2.1 Incorrect spectral modelling of thermal dust 

We now consider a configuration of the sky - Model lb and 
Model 2b - where the thermal dust polarization consists 
in two MBB spectral components (Finkbeiner et al. 1999): 
a cold dust component with temperature Ti = 9.6 K and 
spectral index Pi = 1.5, and a hot dust component with 
temperature T 2 = 16.4 K and spectral index P 2 = 2.6. The 
synchrotron emission still has a uniform power-law distribu¬ 
tion over frequencies, and the CMB either include B-mode 
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Model ID 

Mean 

Mean x^ 

Recovered r 

Experiment 

Model Ob 

37.51 

0.99 

0.05271 ± 0.00595 

COrE+ Light 


41.50 

0.99 

0.05202 ± 0.00585 

COrE+ Extended 

r = 0.05 

29.48 

0.98 

0.05107 ± 0.00575 

COrE 

no foreground 

11.53 

0.96 

0.05132 ± 0.00578 

LiteBIRD 


65.51 

0.99 

0.05145 ± 0.00616 

PIXIE 


13.56 

0.97 

0.05074 ± 0.00572 

EPIC-LC-TES 


15.44 

0.96 

0.05086 ± 0.00583 

EPIC-CS 


17.44 

0.97 

0.05096 ± 0.00572 

EPIC-IM-4K 


39.58 

0.99 

0.05140 ± 0.00578 

PRISM 

Model la 

37.87 

1.00 

0.05224 ± 0.00800 

COrE+ Light 


41.71 

0.99 

0.05132 ± 0.00770 

COrE+ Extended 

r = 0.05 

29.72 

0.99 

0.04877 ± 0.00715 

COrE 

correct foreground 

11.86 

0.99 

0.04997 ± 0.00767 

LiteBIRD 

modelling 

65.67 

1.00 

0.05507 ± 0.00865 

PIXIE 


14.82 

1.06 

0.05411 ± 0.00823 

EPIC-LC-TES 


15.57 

0.97 

0.04506 ± 0.00730 

EPIC-CS 


17.65 

0.98 

0.04989 ± 0.00708 

EPIC-IM-4K 


39.82 

1.00 

0.05027 ± 0.00737 

PRISM 

Model lb 

32.14 

1.07 

0.08929 ± 0.00766 

COrE-l- Light 


35.90 

1.20 

0.09218 ± 0.00624 

COrE-l- Extended 

r = 0.05 

22.32 

1.12 

0.08023 ± 0.01045 

COrE 

incorrect dust 

13.19 

1.10 

0.08428 ± 0.00935 

LiteBIRD 

modelling 

25.86 

1.08 

0.09711 ± 0.00265 

PIXIE 


18.47 

1.32 

0.08113 ± 0.00999 

EPIC-LC-TES 


15.08 

1.08 

0.07911 ± 0.01048 

EPIC-CS 


40.30 

2.88 

0.09434 ± 0.00485 

EPIC-IM-4K 


50.44 

1.58 

0.09446 ± 0.00467 

PRISM 

Model Oa 

37.49 

0.99 

0.00000 ± 0.00005 

COrE-l- Light 


41.49 

0.99 

0.00000 ± 0.00002 

COrE-l- Extended 

r = 0 

29.47 

0.98 

0.00000 ± 0.00004 

COrE 

no foreground 

11.52 

0.96 

0.00000 ± 0.00003 

LiteBIRD 


65.57 

0.99 

0.00000 ± 0.00036 

PIXIE 


13.55 

0.97 

0.00000 ± 0.00003 

EPIC-LC-TES 


15.43 

0.96 

0.00000 ± 0.00022 

EPIC-CS 


17.44 

0.97 

0.00000 ± 0.00001 

EPIC-IM-4K 


39.57 

0.99 

0.00000 ± 0.00001 

PRISM 

Model 2a 

37.83 

1.00 

0.00031 ± 0.00024 

COrE-l- Light 


41.66 

0.99 

0.00016 ± 0.00012 

COrE-l- Extended 

r = 0 

29.72 

0.99 

0.00017 ± 0.00013 

COrE 

correct foreground 

11.85 

0.99 

0.00020 ± 0.00017 

LiteBIRD 

modelling 

65.66 

0.99 

0.00056 ± 0.00042 

PIXIE 


14.88 

1.06 

0.00033 ± 0.00026 

EPIC-LC-TES 


15.57 

0.97 

0.00023 ± 0.00020 

EPIC-CS 


17.60 

0.98 

0.00016 ± 0.00013 

EPIC-IM-4K 


39.78 

0.99 

0.00019 ± 0.00015 

PRISM 

Model 2b 

32.19 

1.07 

0.05229 ± 0.01223 

COrE-l- Light 


35.95 

1.20 

0.06357 ± 0.01332 

COrE-l- Extended 

r = 0 

22.42 

1.12 

0.03453 ± 0.00821 

COrE 

incorrect dust 

13.23 

1.10 

0.01595 ± 0.01249 

LiteBIRD 

modelling 

25.87 

1.08 

0.09246 ± 0.00635 

PIXIE 


18.55 

1.33 

0.02792 ± 0.00936 

EPIC-LC-TES 


15.16 

1.08 

0.03018 ± 0.00699 

EPIC-CS 


40.38 

2.88 

0.07251 ± 0.01265 

EPIC-IM-4K 


50.52 

1.58 

0.06885 ± 0.01068 

PRISM 


Table 3. Maximum likelihood estimate of the tensor-to-scalar ratio r (fourth column) for different experiments for sky simulations with 
baseline r = 0.05 (resp. r = 0) and r = 0.087. Model Ob: no foregrounds, r = 0.05. Model la: correct foregrounds, r = 0.05. Model lb: 
incorrect thermal dust, r = 0.05. Model Oa: no foregrounds, r = 0 (no CMB B-modes). Model 2a: correct foregrounds, r = 0 (no CMB 
B-modes). Model 2b: incorrect thermal dust, r = 0 (no CMB B-modes). In addition, mean values (second column) and normalized 
values (third column) from Bayesian parametric fitting for the separation of the component maps. Nch denotes the number of channels 
for each experiment (twice the number of frequencies because of Q and U polarization maps). 


© 0000 RAS, MNRAS 000, 000-000 



12 M. Remazeilles et al. 


Model ID 

Mean x^ 

Mean x'^l^ch 

Recovered r 

Experiment 

Model Oc 

37.57 

0.99 

0.00104 ± 0.00019 

COrE+ Light 


41.57 

0.99 

0.00103 ± 0.00018 

COrE+ Extended 

r = 0.001 

29.48 

0.98 

0.00104 ± 0.00019 

COrE 

no foreground 

11.56 

0.96 

0.00095 ± 0.00017 

LiteBIRD 


65.56 

0.99 

0.00098 ± 0.00026 

PIXIE 


13.56 

0.97 

0.00097 ± 0.00017 

EPIC-LC-TES 


15.43 

0.96 

0.00114 ± 0.00023 

EPIC-CS 


17.45 

0.97 

0.00100 ± 0.00017 

EPIC-IM-4K 


39.59 

0.99 

0.00106 ± 0.00018 

PRISM 


Table 4. Maximum likelihood estimate of the tensor-to-scalar ratio r (fourth column) for different experiments for a sky simulation with 
baseline r = 0.001 and no foregrounds (Model Oc). In addition, mean values (second column) and normalized values (third column) 
from Bayesian parametric fitting for the separation of the component maps. denotes the number of channels for each experiment 
(twice the number of frequencies because of Q and U polarization maps). 



0.000 0.001 0.002 0.003 0.004 

r (tensor/scolor) 

Figure 2. Recovered posterior distribution P(r) of the tensor- 
to-scalar ratio for Model Oc (theoretical input r = 0.001, 
no foregrounds): COrE (solid yellow), for COrE-f Light (solid 
light-blue), COrE-l- Extended (solid blue), LiteBIRD (dotted 
red), PIXIE (dashed green), EPIC-LC-TES (long-dashed yellow), 
EPIC-CS (long-dashed purple), EPIC-IM-4K (long-dashed or¬ 
ange), and PRISM (dash three-dot black) experiments. 


polarization (r = 0.05 for Model lb) or not {r — 0 for Model 
2b). Therefore, we are considering spectral mismatch be¬ 
tween the simulated thermal dust (two MBB components) 
and the dust model Eq. 7 (one MBB dust component) fitted 
to the data. 

In order to minimize the impact of incorrect dust model¬ 
ling, we exclude from the parametric fitting any frequency 
channel above 360 GHz for all the CMB experiments. This 
is because at higher frequencies the effects of dust temper¬ 
ature variations have a significant effect on the spectrum. 
We therefore restrict the analysis to 10 channels for COrE 
from 45 GHz to 315 GHz, 15 channels for COrE-l- Light from 
60 GHz to 340 GHz, 15 channels for COrE-|- Extended from 
60 GHz to 340 GHz, 6 channels for LiteBIRD from 30 GHz 
to 280 GHz, and 12 channels for PIXIE from 60 GHz to 
360 GHz, 7 channels for EPIC-CS from 30 GHz to 300 GHz, 
7 channels for EPIG-CS from 30 GHz to 340 GHz, 7 channels 
for EPIC-IM-4K from 30 GHz to 340 GHz, 16 channels for 
PRISM from 30 GHz to 320 GHz. 

The results for the recovered distribution of the tensor- 



Frequency [GHz] 

Figure 3. Modified blackbody (MBB) spectral energy distribu¬ 
tions for thermal dust. The two-MBB dust model (solid blue) 
results from the coaddition of a cold dust component (T = 9.6 K, 
/3 = 1.5 - dashed green) and a hot dust component (T = 16.4K, 
/3 = 2.6 — dashed red). The difference between the one-MBB dust 
model (T = 18K, /3 = 1.6 — solid black), adopted in the Bayesian 
fitting, and the two-MBB dust model (solid blue), adopted in 
Model lb and Model 2b sky simulations, is becoming more sig¬ 
nificant at frequencies larger than 300 GHz. 


to-scalar ratio P{r) are shown in the third row of Fig. 1 for 
r = 0 (left-hand panel) and r = 0.05 (right-hand panel). 
Because of incorrect spectral assumptions on the number 
of thermal dust components in the parametric model, the 
recovered tensor-to-scalar ratio is found to be strongly bi¬ 
ased, by more than 3a for all the experiments. The impact is 
even more significant for the experiments having the largest 
sensitivity at high frequencies, i.e. PRISM, EPIG-IM, and 
GOrE-H Extended missions: the recovered tensor-to-scalar 
ratio is in this case biased by more than 5a when incorrectly 
parametrizing the spectral properties of the thermal dust. 
The maximum likelihood estimates of r with Icr errors are 
listed in this case in the third (Model lb) and sixth (Model 
2b) blocks of Table 3, showing strong departure from the 
baseline value of the tensor-to-scalar ratio. The correspond¬ 
ing mean values of the parametric fitting are listed in 
the third and sixth blocks of Table 3: they are significantly 
larger than the baseline values of Model la (second block) 
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Figure 4. Impact of high-frequency channels on the estimation of the tensor-to-scalar ratio r with COrE when incorrectly modelling 
the polarized thermal dust emission. The solid black line is for the frequency range 45-255 GHz and the dashed line is when extending 
the frequency range to 345 GHz. Left panel: Model 2b (r = 0, incorrect dust spectrum). Right panel: Model lb (r = 0.05, incorrect 
dust spectrum). The high-frequency channels of COrE are useful to highlight any failure of the dust model. The instability of the mean 
value of the tensor-to-scalar distribution when including or not high-frequency channels provides a good diagnosis for revealing spurious 
foreground B-mode detection due to incorrect modelling. 


and model 24 (fifth block), for which the foreground model 
was consistent with the data. This indicates a failure in the 
parametrization of the foreground model htted to the data, 
and an unsuccessful separation of CMB and dust B-modes. 

Clearly, the over-prediction of the tensor-to-scalar ra¬ 
tio is due to contamination by spurious dust foreground 
B-modes. The amount of spurious dust B-modes detected 
potentially increases for experiments having more high- 
frequency channels and more sensitivity because spectral 
modelling of thermal dust becomes more and more criti¬ 
cal. Being limited to 255 GHz, LiteBIRD is less impacted 
by incorrect dust modelling, at least when r = 0 (left-hand 
panel of the third row of Fig. 1), because at lower frequen¬ 
cies the coadded spectral distribution of the two MBB dust 
components is still consistent with the single MBB distri¬ 
bution of the model, as illustrated in Fig. 3. However, the 
increased ~ 1-10 value of LiteBIRD (Model 2b in Ta¬ 
ble 3) with respect to the baseline = 0.99 value (Model 
2a in Table 3) still reveals incorrect assumptions on thermal 
dust, in which case the estimate of r can not be considered 
as robust. For PIXIE the bias on r = 0 increases dramat¬ 
ically by more than 14a (Model 2b in Table 3) while the 
-)/ = 1.08 value stays close to the baseline value because 
of the larger noise level in the PIXIE channels compared to 
other experiments. This result shows that low sensitivity in 
a large number of channels may prevent an experiment from 
detecting incorrect foreground modelling, therefore leading 
to a false detection of r. However, in the case of PIXIE, the 
low sensitivity per channel can be compensated by includ¬ 
ing more frequency channels. In summary, experiments with 
more high/low-frequency channels and more sensitivity in 
high/low-frequencies, like PRISM, COrE, COrE-b are better 
designed to control foreground uncertainties and guarantee 
a robust detection of B-modes. 

It is interesting to compare these results with those 
obtained on Planck simulations by Armitage-Caplan et al. 
(2012): modelling a two-component dust simulation with a 
one-component dust model had only a minor impact on 



Figure 5. Top panel: WMAP Low Resolution Polarization Data 
Analysis Mask (Bennett et al. 2013) with /sky = 73 %. Bottom 
panel: Planck-hased polarization mask with /sky = 63%. Blue 
regions are masked. 


Planck, biasing the tensor-to-scalar ratio by less than Icr 
because of relatively lower sensitivity. It is clear that be¬ 
cause next-generation CMB satellite missions (e.g. COrE, 
COrE-b, LiteBIRD, PIXIE, EPIC, PRISM) will have signif¬ 
icantly more sensitivity, they will be much more sensitive to 
any incorrect spectral modelling of thermal dust, as demon¬ 
strated in the third row of Fig. 1. 

The impact on the tensor-to-scalar ratio of incorrect 
dust modelling is more significant for CMB satellites having 
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Figure 6. Impact of Galactic masking on the estimation of the tensor-to-scalar ratio r with COrE when incorrect modelling the thermal 
dust spectrum. The solid black line is for the Planck mask (63 % coverage) and the black dashed line is for the WMAP mask (73 % 
coverage). Left panel: Model 2b (r = 0, incorrect dust model). Right panel: Model lb (r = 0.05, incorrect dust model). The Planck mask, 
designed for masking the thermal dust intensity distribution over sky, reduces the amount of bias in the estimated tensor-to-scalar ratio. 


high-frequency observations. This is due to the increasing 
significance of the polarized dust emission at high frequency, 
where spectral parametrization becomes critical for this fore¬ 
ground. From the results shown on the third row of Fig. 1, 
it is clear that high-frequency channel observations (350 to 
850 GHz) are particularly useful for CMB satellite missions 
to highlight any failure in the modelling of the polarized 
thermal dust emission. 

Taking the COrE satellite as an example, Fig. 4 shows 
the effect of frequency coverage on the bias on r when in¬ 
correctly parametrizing the thermal dust foreground. The 
amount of bias in the recovered tensor-to-scalar ratio with 
respect to the original value r = 0 (left-hand panel) or 
r = 0.05 (right-hand panel) increases when including in the 
parametric fitting frequencies larger than 255 GHz. 

The check of instability of the mean value of the tensor- 
to-scalar distribution P(r) when including or not higher fre¬ 
quency channels in the Bayesian analysis, as shown in Fig. 4, 
provides a robust test of spurious foreground B-mode detec¬ 
tion and incorrect foreground modelling. 

The amount of bias of the recovered tensor-to-scalar ra¬ 
tio due to incorrect foreground modelling is also dependent 
on Galactic masking. To illustrate this, we use two different 
Galactic masks (Fig. 5): the WMAP Low Resolution Polar¬ 
ization Data Analysis Mask® (Bennett et al. 2013), relying 
mainly on synchrotron intensity thresholds over the sky from 
the WMAP K band polarized intensity map, and a more 
conservative polarization mask, based on both synchrotron 
and thermal dust intensity thresholds from the Planck po¬ 
larized intensity maps at 30 GHz and 353 GHz (Planck Gol- 
laboration et al. 2015a). The synchrotron intensity thresh¬ 
old is determined as follows. We simulate GMB Q and U 
polarization maps at 70 GHz and 5 degree beam resolution 
and we compute the RMS of the polarized intensity of the 
GMB, — 0.4 ^K, which is taken as a reference. 

We smooth the Planck 30 GHz polarization intensity map at 


® http;//lambda.gsfc.nasa.gov/product/map/dr5/masks_get. 
cfm 


5 degree resolution and extrapolate it to 70 GHz through a 
power-law with spectral index /3 = —3. We compute its RMS 
and mask any pixel where the RMS of the Planck 30 GHz 
polarization intensity extrapolated to 70 GHz is 2.5 larger or 
more than the reference RMS of the GMB polarization inten¬ 
sity. The dust intensity threshold is determined in the exact 
same way by using the Planck 353 GHz polarization map 
and extrapolating it throuh a power-law with spectral index 
P = 1.6. The Planck-ha.seA conservative mask is defined by 
the union of the synchrotron mask and the dust mask. The 
fraction of the observed sky is 63 % for the Planck-haseA 
conservative mask and 73 % for the WMAP mask. 

Figure 6 shows the effect of masking for the COrE ex¬ 
periment when incorrectly modelling the thermal dust po¬ 
larization. The bias in the recovered cosmological signal r is 
reduced when using a mask that is designed from thermal 
dust intensity thresholds, such as the Planck-hased mask. 
Masking the Galactic plane 10 % more conservatively leads 
to a significant change in the bias of the recovered tensor-to- 
scalar ratio. Again, a check of instability of the mean value of 
the tensor-to-scalar distribution when changing the Galac¬ 
tic mask used in the component separation process provides 
a robust test of spurious foreground B-mode detection and 
incorrect foreground modelling. 

This highlights a major advantage of the parametric 
pixel-fitting method employed here: in practice, the map 
can be used to update the model for pixels that do not pro¬ 
vide a good fit to the data, or, alternatively, the mask can 
be increased to ignore those pixels thereby improving the 
overall result (Eriksen et al. 2008). 


5.8.2 Incorrect spectral modelling of synchrotron 

While polarization of the Galactic synchrotron emission is 
expected to scale with frequency approximately according 
to a power-law spectral distribution, with Ps ~ —3 
(Planck Gollaboration et al. 2015d), the exact value of the 
synchrotron spectral index ps {p, v) is certainly non-uniform 
over the sky (e.g., Platania et al. 1998; Miville-Deschenes 
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Model ID 

Mean 

Mean 

Recovered r 

Experiment 

Model Ic 

37.88 

1.00 

0.05343 ± 0.00830 

COrE+ Light 


41.73 

0.99 

0.05212 ± 0.00777 

COrE+ Extended 

r = 0.05 

29.78 

0.99 

0.05049 ± 0.00741 

COrE 

non-uniform 

11.88 

0.99 

0.05094 ± 0.00799 

LiteBIRD 

synchrotron 

65.69 

1.00 

0.05620 ± 0.00866 

PIXIE 

spectral index 

14.96 

1.07 

0.05507 ± 0.00827 

EPIC-LC-TES 


15.65 

0.98 

0.04790 ± 0.00766 

EPIC-CS 


18.01 

1.00 

0.05110 ± 0.00726 

EPIC-IM-4K 


40.14 

1.00 

0.05075 ± 0.00745 

PRISM 

Model Id 

38.15 

1.00 

0.06756 ± 0.01027 

COrE+ Light 


42.29 

1.01 

0.06390 ± 0.00946 

COrE+ Extended 

r = 0.05 

30.30 

1.01 

0.06074 ± 0.00920 

COrE 

incorrect 

12.14 

1.01 

0.07988 ± 0.01027 

LiteBIRD 

synchrotron 

66.47 

1.01 

0.07122 ± 0.01027 

PIXIE 

curvature C=+0.3 

15.21 

1.09 

0.07769 ± 0.01029 

EPIC-LC-TES 


15.91 

0.99 

0.06558 ± 0.01004 

EPIC-CS 


23.43 

1.30 

0.06205 ± 0.00906 

EPIC-IM-4K 


44.82 

1.12 

0.06386 ± 0.00925 

PRISM 

Model le 

38.07 

1.00 

0.05784 ± 0.00888 

COrE-l- Light 


42.21 

1.01 

0.05500 ± 0.00817 

COrE-l- Extended 

r = 0.05 

30.53 

1.02 

0.05318 ± 0.00780 

COrE 

incorrect 

12.17 

1.01 

0.06356 ± 0.00972 

LiteBIRD 

synchrotron 

65.84 

1.00 

0.06119 ± 0.00939 

PIXIE 

curvature C=-0.3 

15.65 

1.12 

0.07379 ± 0.01030 

EPIC-LC-TES 


16.04 

1.00 

0.05817 ± 0.00899 

EPIC-CS 


21.27 

1.18 

0.05267 ± 0.00750 

EPIC-IM-4K 


43.02 

1.08 

0.05263 ± 0.00773 

PRISM 

Model 2c 

37.83 

1.00 

0.00067 ± 0.00028 

COrE-l- Light 


41.68 

0.99 

0.00018 ± 0.00014 

COrE-l- Extended 

r = 0 

29.79 

0.99 

0.00021 ± 0.00016 

COrE 

non-uniform 

11.88 

0.99 

0.00076 ± 0.00035 

LiteBIRD 

synchrotron 

65.67 

1.00 

0.00086 ± 0.00057 

PIXIE 

spectral index 

15.01 

1.07 

0.00059 ± 0.00038 

EPIC-LC-TES 


15.65 

0.98 

0.00047 ± 0.00030 

EPIC-CS 


18.01 

1.00 

0.00021 ± 0.00017 

EPIC-IM-4K 


40.13 

1.00 

0.00023 ± 0.00018 

PRISM 

Model 2d 

38.11 

1.00 

0.01252 ± 0.00203 

COrE-l- Light 


42.25 

1.01 

0.01019 ± 0.00151 

COrE-l- Extended 

r = 0 

30.27 

1.01 

0.01024 ± 0.00162 

COrE 

incorrect 

12.19 

1.02 

0.02608 ± 0.00428 

LiteBIRD 

synchrotron 

66.46 

1.01 

0.01313 ± 0.00236 

PIXIE 

curvature C=+0.3 

15.24 

1.09 

0.01830 ± 0.00283 

EPIC-LC-TES 


15.93 

1.00 

0.01647 ± 0.00269 

EPIC-CS 


23.45 

1.30 

0.01135 ± 0.00148 

EPIC-IM-4K 


44.83 

1.12 

0.01211 ± 0.00156 

PRISM 

Model 2e 

38.03 

1.00 

0.00499 ± 0.00107 

COrE-l- Light 


42.15 

1.00 

0.00338 ± 0.00062 

COrE-l- Extended 

r = 0 

30.48 

1.02 

0.00495 ± 0.00131 

COrE 

incorrect 

12.21 

1.02 

0.01149 ± 0.00260 

LiteBIRD 

synchrotron 

65.82 

1.00 

0.00538 ± 0.00136 

PIXIE 

curvature C=-0.3 

15.72 

1.12 

0.01829 ± 0.00258 

EPIC-LC-TES 


16.03 

1.00 

0.01080 ± 0.00224 

EPIC-CS 


21.22 

1.18 

0.00322 ± 0.00042 

EPIC-IM-4K 


42.98 

1.07 

0.00324 ± 0.00093 

PRISM 


Table 5. Maximum likelihood estimate of the tensor-to-scalar ratio r (fourth column) for different experiments for sky simulations with 
baseline r = 0.05 (resp. r = 0) and r = 0.087. Model Ic: neglected variability of synchrotron spectral index, r = 0.05. Model Id: neglected 
positive synchrotron curvature C = +0.3, r = 0.05. Model le: neglected negative synchrotron curvature C = —0.3, r = 0.05. Model 
2c: neglected variability of synchrotron spectral index, r = 0 (no CMB B-modes). Model 2d: neglected positive synchrotron curvature 
C = +0.3, r = 0 (no CMB B-modes). Model 2e: neglected negative synchrotron curvature C = —0.3, r = 0 (no CMB B-modes). In 
addition, mean values (second column) and normalized x^ values (third column) from Bayesian parametric fitting for the separation 
of the component maps. Nch denotes the number of channels for each experiment (twice the number of frequencies because of Q and U 
polarization maps). 
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Figure 7. Recovered posterior distribution P{r) of the tensor-to-scalar ratio and impact of incorrect synchrotron modelling. The 
theoretical input tensor-to-scalar value (vertical solid black line) is r = 0 in the left-hand panels and r = 0.05 in the right-hand panels. 
Top panels-, neglected variability of the synchrotron spectral index {left-. Model 2c, right: Model Ic). Middle panels: neglected synchrotron 
curvature C = +0.3 {left: Model 2d, right: Model Id). Bottom panels: neglected synchrotron curvature C = —0.3 {left: Model 2e, right: 
Model le). Recovered tensor-to-scalar distributions: COrE (solid yellow), COrE+ Light (solid light-blue), COrE+ Extended (solid blue), 
LiteBIRD (dotted red), PIXIE (dashed green), EPIC-LC-TES (long-dashed yellow), EPIC-CS (long-dashed purple), EPIC-IM-4K (long- 
dashed orange), PRISM (dash three-dot black). 


et al. 2008; Irfan et al. 2015) and may also vary with fre¬ 
quency (curvature) (Kogut et al. 2007; de Oliveira-Costa 
et al. 2008; Kogut 2012). These uncertainties on the spec¬ 
tral properties of the synchrotron may result in an incorrect 
fitting to the data at low frequency, therefore an incorrect 
foreground subtraction, and subsequently a biased estima¬ 
tion of the tensor-to-scalar ratio. 

In the top panels of Fig. 7 (Models Ic and 2c), we con¬ 


sider a configuration where the sky data include synchrotron 
polarization with a non-uniform spectral index over the sky 
(Miville-Deschenes et al. 2008) but we parametrize the syn¬ 
chrotron in the fitting model by a uniform /3s = — 3 ± 0.1 
spectral index (Gaussian prior). The complete set of fre¬ 
quency channels is used in the parametric fitting for all 
CMB experiments. Making the wrong assumption of a uni¬ 
formly distributed synchrotron spectral index for r = 0.05 
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Figure 8. Recovered posterior distribution P{r) of the tensor-to-scalar ratio and impact of neglecting the spinning dust polarization. The 
theoretical input tensor-to-scalar value (vertical solid black line) is r = 0 in the left-hand panel (Model 21) and r = 0.05 in the right-hand 
panel (Model If). Recovered tensor-to-scalar distributions: COrE (solid yellow), COrE-|- Light (solid light-blue), COrE-|- Extended (solid 
blue), LiteBIRD (dotted red), PIXIE (dashed green), EPIC-LC-TES (long-dashed yellow), EPIC-CS (long-dashed purple), EPIC-IM-4K 
(long-dashed orange), PRISM (dash three-dot black). 


(top right-hand panel of Fig. 7) has a negligible impact on 
the recovered tensor-to-scalar ratio for all CMB satellites, 
the corresponding estimates of r listed in the first block of 
Table 5 show consistent values with the baseline estimates of 
Model la (Table 3) where the foreground model was match¬ 
ing the data. Therefore, at the level of sensitivity of the CMB 
satellites considered in this work, the non-uniform distribu¬ 
tion over the sky of the synchrotron spectral index (Miville- 
Deschenes et al. 2008) is successfully fitted by the Gaussian 
distribution /3s = —3±0.1 on /sky = 73% of the sky (Fig. 5). 

When r = 0 (top left-hand panel of Fig. 7), COrE-b 
Light and LiteBIRD experiments recover a non-zero value of 
the tensor-to-scalar ratio. However, the mean of the para¬ 
metric fitting (Model 2c in Table 5) for these experiments 
show consistent values with the baseline Model 2a (correct 
foreground modelling), therefore indicating a good fit to the 
data in the frequency range considered. In this case, the 
bias on the tensor-to-scalar ratio is rather attributed to the 
lack of low-frequency channels for COrE-b Light and Lite¬ 
BIRD, for which iz ^ 60 GHz. The lack of low-frequency ob¬ 
servations aX V < QQ GHz for GMB satellites prevents them 
to fit for a distribution of the synchrotron spectral index 
that is more complex than a simple constant spectral index. 
As shown in top left-hand panel of Fig. 7 and the fourth 
block of Table 5, because of the absence of low-frequency 
channels, LiteBIRD and COrE-b Light missions may fail in 
recovering the correct tensor-to-scalar ratio if the spectral 
distribution of the polarized synchrotron is non-trivial at 
low frequency. Sensitive low frequency surveys (e.g. C-BASS 
at 5 GHz; Irfan et al. 2015) and also the lowest frequency 
channels of WMAP/Plancfc may also help. For example, for 
the case of LiteBIRD, including a C-BASS data point with 
O.lmKbeam"^ noise reduces the final CMB uncertainty by 
~ 50 %; of course the exact value is very sensitive on the 
foreground model adopted. 

In the middle panels (resp. bottom panels) of Fig. 7, 
we show the impact on the estimated tensor-to-scalar ratio 
of neglecting in the foreground model a C = -bO.3 (resp. 


C = —0.3) curvature of the synchrotron spectral index.® 
Neglecting a positive synchrotron curvature significantly bi¬ 
ases the recovered tensor-to-scalar ratio by Icr at least for all 
CMB experiments (middle right-hand panel of Fig. 7). Note 
that the additional bias found for LiteBIRD compared to 
others is not due to incorrect fitting but again to the lack of 
low-frequency channels where synchrotron modelling is crit¬ 
ical. Comparing the middle and bottom right-hand panels 
of Fig. 7 shows that neglecting a positive curvature of the 
synchrotron spectral index has more impact on the tensor- 
to-scalar estimation that omitting a negative synchrotron 
curvature. This can be interpreted by the fact that flattening 
the synchrotron distribution at increasing frequency makes 
the synchrotron less “orthogonal” to the CMB spectral dis¬ 
tribution and makes the separation of CMB and synchrotron 
more challenging. 

5.2.3 Neglecting spinning dust polarization 

The number of polarized foregrounds in the sky is not per¬ 
fectly known. Apart from thermal dust and synchrotron 
emissions, extra polarized components such as free-free, 
magnetic dust (MD), and spinning dust emissions should 
not be omitted. Free-free emission is probably not a major 
component since it is intrinsically unpolarized (Rybicki & 
Lightman 1986; Keating et al. 1998; Macellari et al. 2011) 
while magneto-dipole radiation from magnetized dust grains 
could be significantly polarized, up to « 35 % according to 
some theoretical models (Draine & Lazarian 1999; Draine & 
Hensley 2013). The latter could potentially be a major CMB 
polarization foreground at ~ 70-150 GHz where the highest 
sensitivity channels are located (see Table 1). There has yet 
to be a definitive detection of MD although there are several 
indications that it may be contributing significantly to the 
total emission at frequencies « 50-300 GHz (Draine & Hens¬ 
ley 2012; Planck Gollaboration et al. 2015f) and thus it will 

® We adopt the convention that positive curvature corresponds to 
a flattening of the spectral distribution of the synchrotron emis¬ 
sion, as given in Equation 3. 
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Model ID 

Mean 

Mean x^l^^ch 

Recovered r 

Experiment 

Model If 

37.89 

1.00 

0.05307 ± 0.00824 

COrE-l- Light 


41.77 

0.99 

0.05159 ± 0.00767 

COrE+ Extended 

r = 0.05 

30.63 

1.02 

0.05158 ± 0.00767 

COrE 

omitted 

11.93 

0.99 

0.05178 ± 0.00791 

LiteBIRD 

spinning dust 

65.72 

1.00 

0.05606 ± 0.00861 

PIXIE 


16.81 

1.20 

0.06956 ± 0.01015 

EPIC-LC-TES 


16.05 

1.00 

0.05208 ± 0.00826 

EPIC-CS 


24.76 

1.38 

0.05165 ± 0.00737 

EPIC-IM-4K 


51.18 

1.28 

0.05191 ± 0.00762 

PRISM 

Model 2f 

37.85 

1.00 

0.00044 ± 0.00025 

COrE-l- Light 


41.73 

0.99 

0.00017 ± 0.00013 

COrE+ Extended 

r = 0 

30.67 

1.02 

0.00051 ± 0.00036 

COrE 

omitted 

11.93 

0.99 

0.00026 ± 0.00022 

LiteBIRD 

spinning dust 

65.70 

1.00 

0.00068 ± 0.00045 

PIXIE 


16.95 

1.21 

0.00573 ± 0.00118 

EPIC-LC-TES 


16.08 

1.01 

0.00155 ± 0.00051 

EPIC-CS 


24.88 

1.38 

0.00026 ± 0.00021 

EPIC-IM-4K 


51.29 

1.28 

0.00032 ± 0.00024 

PRISM 


Table 6. Maximum likelihood estimate of the tensor-to-scalar ratio r (fourth column) for different experiments for sky simulations with 
baseline r = 0.05 (resp. r = 0) and r = 0.087. Model If: omitted spinning dust, r = 0.05. Model 2f: omitted spinning dust, r = 0 (no 
CMB B-modes). In addition, mean values (second column) and normalized values (third column) from Bayesian parametric fitting 
for the separation of the component maps. denotes the number of channels for each experiment (twice the number of frequencies 
because of Q and U polarization maps). 


need to be considered carefully when applying component 
separation. 

The Galactic emission from spinning dust grains on the 
other hand is known to be significant with several clear de¬ 
tections at least along some sight-lines (Planck Collabora¬ 
tion et al. 2011b, 2015a) and theoretical models suggest it 
should be slightly polarized at a fraction of 1 % or less (see 
Sect. 3). Given the expected sensitivity of future CMB satel¬ 
lite missions, the small B-mode polarization of the spinning 
dust may play a non-negligible role when estimating the 
tensor-to-scalar ratio. 

In Fig. 8 we show the impact on the tensor-to-scalar 
ratio of neglecting in the foreground model the 1% polar¬ 
ization of spinning dust. The sky simulations Model If and 
Model 2f differ from the baseline simulations Model la and 
Model 2a just by the insertion of spinning dust polarization 
in the data. However, when omitting the 1 % polarization 
of the spinning dust in the foreground model, the recov¬ 
ered tensor-to-scalar ratio is clearly impacted as it is shown 
by comparing the panels of Fig. 8 (Models If and 2f) with 
the second row of panels in Fig. 1 (Models la and 2a). In 
particular, neglecting spinning dust polarization results in 
a non-negligible bias on the recovered tensor-to-scalar ratio 
for the EPIC experiment (also see Table 6), with a large 
(Table 6) compared to the baseline value, thus indicating 
imperfect component separation. Typically, neglecting the 
spinning dust has more impact on the CMB satellites hav¬ 
ing larger sensitivity at low frequency (« 30-45 GHz), where 
the spinning dust polarization intensity is more significant. 


6 CONCLUSIONS 

The Bayesian framework for B-mode component separation 
allows the propagation of errors in foreground modelling to¬ 


wards the estimation of the tensor-to-scalar ratio. It is also 
useful to re-adjust, a posteriori, the model of Galactic fore¬ 
ground polarization being fitted by using feedback from the 
map of the goodness-of-fit statistics. 

We have compared the overall sensitivity to B-modes of 
different CMB satellite concepts both in the absence and in 
the presence of foreground contamination. This is the first 
time that these proposed experiments have been compared 
on an equal basis. 

Due to unprecedented sensitivity of next-generation 
CMB satellites, any incorrect spectral assumptions on the 
Galactic foregrounds may result in a significantly biased es¬ 
timation of the tensor-to-scalar ratio after foreground clean¬ 
ing. Modelling two MBB thermal dust components as a sin¬ 
gle MBB thermal dust component strongly biases the es¬ 
timate of the tensor-to-scalar ratio by more than 5a for 
the most sensitive experiments like PRISM, EPIC-IM, and 
GOrE-k extended. Neglecting the synchrotron curvature, 
C = 0.3, biases all CMB satellites by la to 2a. Omitting 
the 1 % polarization of spinning dust in the parametrization 
of Galactic foregrounds makes a non-negligible bias on the 
recovered tensor-to-scalar ratio. 

In this work we have used two distinct criteria to eval¬ 
uate the accuracy of the CMB B-mode estimation: the x^ 
statistics of the parametric fitting and the bias on the re¬ 
covered tensor-to-scalar ratio. Taken together these crite¬ 
ria indicate whether the incorrect estimation of the tensor- 
to-scalar ratio results from incorrect fitting of the Galactic 
foregrounds or from a lack of frequency channel observa¬ 
tions. In particular, the lack of high/low-frequency channels 
(e.g. LiteBIRD) may prevent from fitting for non-trivial fore¬ 
grounds. Most important, we have shown how the lack of 
high/low-frequency channels and the low sensitivity in indi¬ 
vidual channels for can lead to a low x^ value from para¬ 
metric fitting with still a biased estimate of r, therefore pre- 
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venting such setups from detecting the incorrect foreground 
modelling and leading to false detection of r. A positive cur¬ 
vature of the synchrotron spectral index flattens the spec¬ 
trum of the synchrotron, which over a limited frequency 
range makes it less distinct from the flat CMB spectrum. 
While the fit of the global sky emission can be accurate 
over the restricted frequency range resulting in a low 
yet the CMB and the curved synchrotron are not success¬ 
fully splitted, therefore leading to an extra bias on r. Note 
that an extended setup of the LiteBIRD mission including 
extra frequency channels and broader frequency coverage 
will be proposed to JAXA (Masashi Hazumi and Tomotake 
Matsumura, priv. comm.). Such an extension of the original 
LiteBIRD concept goes in the right direction of the conclu¬ 
sions of the paper. 

We have limited the analysis to low multipoles {£ < 12) 
where leasing B-modes are not expected to dominate the 
uncertainty on the tensor-to-scalar ratio. However, if r is 
low and r is low (in particular, latest results from Planck 
Collaboration et al. (2015e) indicate a lower value r = 0.066) 
then the error due to leasing B-modes, more than cosmic 
variance, should be considered in addition to foregrounds 
for a complete comparison of the satellite concepts (Smith 
et al. 2012). While a large panel of frequency channels is 
helpful for CMB satellites to clean non-trivial foreground 
B-modes, high-resolution channels should also be useful to 
improve the ability of the experiment to correct for leasing 
B-modes. 

Given the expected sensitivity of next-generation CMB 
satellite missions, it appears to be fundamental to accurately 
characterize the spectral properties of the Galactic fore¬ 
grounds in order to detect primordial CMB B-modes and to 
correctly measure the tensor-to-scalar ratio. In this respect, 
CMB satellite missions with both low-frequency (< 60 GHz) 
and high-frequency (> 300 GHz) channels will have enough 
degrees of freedom to characterize the foregrounds. Sensi¬ 
tive foreground-dedicated experiments at low (< 15 GHz) 
frequencies such as S-PASS (Garretti 2011), G-BASS (Irfan 
et al. 2015) and QUIJOTE (Rubino-Martln et al. 2012b), 
and balloon experiments at high (> 300 GHz) frequencies 
such as PILOT (Misawa et al. 2014) will also be very help¬ 
ful. 
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